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Abstract 

The differential equation in the external invariant satisfied by the master 
integral of the general massive 2-loop 4-denominator self-mass diagram is exploited 
and the expansion of the master integral at = is obtained analytically. The 
system composed by this differential equation with those of the master integrals 
related to the general massive 2-loop sunrise diagram is numerically solved by the 
Runge-Kutta method in the complex p"^ plane. A numerical method to obtain 
results for values of at and close to thresholds and pseudo-thresholds is discussed 
in details. 
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1 Introduction. 



Precise measurements of particles properties require that the corresponding theoretical 
calculations have to include up to (at least) two-loop radiative corrections. In this paper 
we investigate a fast and flexible numerical method for their accurate evaluation. 

A commonly used procedure in modern radiative correction calculations is to express 
the result as a combination of a limited number of Master Integrals (Mi's), using the 
integration by parts identities pQ. In this framework, it is necessary to obtain a precise 
determination of the Mi's, even if the analytical values cannot be obtained due to the 
large number of different scales occurring in each of the Mi's (internal masses and external 
momenta or Mandelstam variables), as it happens in electroweak theory. 

The general massive 2-loop self-mass diagram has four Mi's related to the sunrise (3- 
denominator) diagram, only one independent MI related to the 4-denominator diagram, 
and again only one new MI related to the 5-denominator diagram |2l El IH E] • 

Several procedures for a precise numerical evaluation of the Mi's have been and still 
are investigated, such as multiple expansions jB] , numerical integration jHl 13 El 13 UHl UH 
1121 1131 mj 5 or difference equations [1^] . 

Another method based on differential equations fHl 1^ was proposed in ^H] , where it 
was shown how to use Runge-Kutta method ^Hl to solve numerically the system of linear 
differential equations in the external invariant satisfied by sunrise master integrals. The 




Figure 1: The general massive 2-loop 4-denominator self- mass diagram. 

method was extended to 2-loop 4-denominator and 5-denominator cases in [2015 where it 
was also suggested how to evaluate numerically the Mi's nearby thresholds and pseudo- 
thresholds. We give in this paper an accurate implementation of that approach for the 
sunrise and the 4-denominator Mi's. 

In Section 2, the (n — 4) expansions of the Mi's are constructed, in Section 3 initial 
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conditions for the differential equations for the Master Integrals (or Master Differential 
Equations, MDE's) are discussed. Some results of the program, showing characteristic 
behaviour of the studied Mi's are presented in Section 4. In Section 5, a method for the 
numerical evaluation of the Mi's near thresholds and pseudo-thresholds is discussed in 
detail, while Section 6 is devoted to comparisons with existing results. 

Finally, in Section 7, our conclusions on the application of the method to present and 
further work are presented. 



2 The MDE and the expansion in (n — 4) of the MI 

We use here the following definition of the MI related to the general massive 2-loop 4- 
denominator self-mass diagram, shown in Figd 



{kf + mj) [{p - kiY + m|] (A;| + ml) [{p - ki - ^2)^ + ml] 



(1) 



where integration is performed in n— dimensional Euclidean space. ^ 

Wherever necessary to avoid ambiguities, the usual imaginary displacements in the 
masses mf — >■ mf — ie, where e is an infinitesimal positive number, are understood. The 
arbitrary mass scale /i accounts for the continuous value of the dimensions n. In numerical 
calculations, we choose /i = mi -|- m2 + m^, one of the natural scales of the problem, 
corresponding to the 3-body threshold, while for simplicity in all analytic formulae we 
put ji = 1. To recover results for arbitrary /x, one has to substitute rrii/fi and 

p2 ^ p'^/fi^. 

The master equation reads * 

d_ 

Qp2 

" r>2/ 2 2 2 \ i 2 2 2 2 2\ 

-R [-p ,m^,mjFi{n,m^,m2,m^,m^,p J 



T->2/ 2 2 2\ 2 ^ / 2 2 2 2 2\ 

R [-p .m^.m^) p —Fi{n,m^,m2,m^,m4^,p 
n - 4 



2 

/ o\ [/ 2 I 2\ 2 I / 2 2\2l ^^ / 2 2 2 2 2 

— {n — 6) [m^ + m^)p + [m^ — m^) J:'4[n,m^,m2,m^,m4,p 

, /o 2 2 I 2\ 2 T-. / 2 2 2 2\ 

+ [op — m^ + m^)mil'i{n,m^,m2,m^,p ) 
3n 



+ {p - m^ + m^ 



771 / 2 2 2 2 

-Fo{n,m^,m2,m^,p 



'^F4{n, ml,m2, m§, m1,p'^) of the present paper corresponds to C'^{n)G{n, mf^ml, to|, ml,p'^) of 
■'•Note the change in sign in the 3rd Hne of Eq.(|2| in comparison to [5], due to a misprint there. 



2 



I 2 Tn f 2 2 2 2\i 2 f 2 2 2 2\ 

- ^(n - 2)V{n, ml, m^, ml) 



(2) 



where 



r>2f 2 2 2\ 4, 4, 4io22|o22 o22 

/t (— p ,m-^,m^) = p +mi+m^ + 2mip + Zm^p — Zm^m^ 

= [p"^ + {mi + m4,y][p'^ + {mi - m^)^] . (3) 

The Fj{n,m1,m2,ml,p'^),j = 0,1,2,3 are the massive 2-loop sunrise self-mass master 
amphtudes, discussed in jlj and numerically calculated in |18|, 

Fj{n,mlml,mlp^) = ^^2^)^.2)2 j dJ'ki j (rk2 

^ J =0,1,2,3 (4) 



{kl + mf)"i(j)(fc| + m|)°2(j)((p ~ki- k2Y + m|)°3(j) 



where: for j = 0, i = 1,2,3, ai(0) = 1; for j,i = 1,2,3, ai{j) = 1, if i 7^ j; ai{j) = 2, 
for i = j. The function V{n,mi,ml,m'^), finally, is the massive 2-loop sunrise vacuum 
amplitude defined in Q (see also j^H I 



,8-2n 



V{n,ml,ml,ml) = / d^ki J ^ ^ ^^^^ ^ ^^^^ ^ . (5) 

In the following we also use the massive 1-loop self-mass 

S{n,m\,ml,p^) = ( (Tk J— ^r—^ — —— 3-, (6) 



(27r)"-2 J (P + ml) [{p - + m|] ' 

with minor changes in the notation with respect to [T7j . 
The expansion in {n — 4) of the solution of Eq. (j21) reads 



2 2 2 2 2\ ri2t \) ^ 17(^2)/ 2 2 2 2 2\ 

i,m2,m3,m4,p ) = C {n)i^——^Fl ' {m-^,m2,m^,m^,p ) 

(^l)('^2 2 2 2 2\ I 17(0)/ 2 2 2 2 2\ 



+ (^731)^4 ''{mi,mimimip')+F^'"{mi,mimimip') + 0{n-4)^ , (7) 



where the coefficient 



C(„) = (2v/?)'*-"'r(3-0 



is not expanded to simplify analytical results. Similar Laurent-expansions in {n — 4) 
of Fj{n,ml,ml,ml,p'^), (for j = 0,1,2,3) and l^(n, , m^, m|) were presented in |4j, 
where explicit analytic expressions for the expansion coefficients F- '^\m\,m%m'^,p^), 
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Fj ^\m1,ml,ml,p'^),V'- ^^(m^, m^, mg), V^^ -'^^(m^, mg) and V^(°)(m^, mj, mg) are also 
given, while the Fj^\m'l, m^, m^^p^) (for j = 0, 1, 2, 3) can be found numerically following 



the algorithm outlined in |T8]. By substituting the expansions of all the amplitudes in 
the MDE Eq.(j2I) and using the results of the first coefficients are found to be 



.(-2), 



,2 2\ 



1 



2 „2\ 



mi,m2,mg,m4,p 



lb 



(9) 



where S^^\m1, ml,p'^) is the finite part, at n = 4, of the expansion of the MI from Eq.| 



S{n,mj^,m^,p ) = C{n) 



2{n-A) 



+ S^'^\ml,mlp^) + 0{n-4) 



(10) 



Its analytic value is 



log(mim4) 



4p2 

where R{—p'^,'m1,ml) 



R{-p\ ml ml) \og{u{p\ ml ml)) + (m? - ml) log 



mi 

7714 



B? {—p"^ , ml, m\) and 



/ 2 2 2 

u(p , m^^, m4 



^Jp^ + {mi + 7714) 2 — ^Jp^ + {mi — m^Y" 



p^ + (mi + m4)2 + Jp^ + (mi — m4) 



(11) 



The function F^\m\, ml ml m|,p^) is not known analytically, but it satisfies the differ- 
ential equation 



d 



r>2/ 2 2 2n 2 17(0)/ 2 2 2 2 2n / 2 , 2n 217(0)/ 2 2 2 2' 

R {-p ,m^, m^) p ^-^4 ("^l,"^2 5"^35"^4'P ) = + '^a)'P -^4 ("^l, "^2? "^35 "^4'P 



/ 2 2\2 f 17(0)/ 2 2 2 2 2n 1 o(0) / 2 2 2\ 

■(m^ -m^) '{m^,m2,m^,m^,p ) - -S'^ >{mi,m^,p ) 



-{ml - ml) 



n T^iO) / 2 2 2 2\ I 2 77(UI/ 2 2 2 



2 17(0), 



,2 ^2n 



I 217(0)/ 2 2 2 2\ I 277(0)/ 2 2 2 2\ t/(0)/ 2 2 
+m2-r2 1^1? ""^2? ""^S^P J + ^3-^3 {''^11^21 '^a,iP ) ~ ^ {^2} '^3,1 



+ - 



mj 



4 



/ 2 I 2\ I 2 

(m2 + mg) + -m4 



mi 



+ 



m4 
16" 



[mi — mi) 
8 



3(m2 + mg) — m; 

1 



-mi log(mi) + m2 log(m2) + mg log(mg) - -m4 log(m4) 
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ori(O)/' 2 2 2 2\ I o 2riiui/ 2 2 2 2\ , 2 7-^^U)/ 2 2 2 2\ 



.2 r{0) 



.2ir(0), 



9 7— rfO) / 9 9 9 9\ T / 9 9 9\ 15 o ^ / 9 

o4 io 



1 



1 



2 „2^ 



:i2) 



1 1 

.64 ^ 4' 

Although the coefficient S^^\m1, ml, p'^) is given analytically above, we report here also 
the differential equation, as it is used in the numerical program, for the related quantity 

S(0)(m2,ml,p2) = ^(0)(m2,mip2) _5(o)(^2^^2^0), 

d ^, 



,2 ^2\ 



P^^S^^\ml,ml,p- 



dp"^ 



-S'^'\mlmiy) 



+ 



p 



1 2\ 

m^,p ) 



\p^ + (mi — m4)2] 2 



(mf 



m\m^ I m 

log 



m|) 



ml 



mim4 



/ 2 2xlog(-T 



(13) 



3 Initial conditions (expansion at 7)^ ~ ). 



To start the Runge-Kutta advancing solution method we choose as initial point the special 
point p^ = 0, which allows the analytic calculation. However, as this is one of the points 
where the coefficient multiplying the derivative of ^4"^ vanishes (Eq. lfT^ ). it is necessary 
to know also the second term in the expansion at ~ 0. 

The general massive 1-loop self-mass expansion at p"^ ^ 0, was presented in but 
we report here the explicit formulae to uniform the notation 

S^'\mlmlp') = S!,'\mlml) +p's['\mlml) +0{{py) , (14) 



where 

sS'\mlml)^S^'\mlmlO) 



[mf — m\) 

1 1 
m\m\ 



log(mi) - 



m; 



m^ 



■ log(m4) - 1 



(mf 



m 



2^3 



m^ + m4 



m2 



2(mf 



m 



2\2 



(15) 



As = is a regular point, we define the expansion at ~ of the MI from Eq.(^ as 

77/ 2222 2\ T? ( 222 2\ 

F^{n,m^,m2,m^,m^,p ) = F^^o{n,m^,m2,m^,m^) 

+ F4^,{n,mlmlmlml) p^ + 0{{p''f) (16) 
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The value at = is easily found from Eq.(ITI) and reads 



F4^fi{n,m'l,ml,ml,m1) = F4^{n, ml, ml, 1711,1711,0) 

V{n, ml, ml, ml) — V{n, m\, ml, m^) 



1 1 
m\ — m\ 



It can be in turn expanded in (n — 4) in the usual way 



with 



F4,o{n,ml,ml,ml,ml) = C'\n)^j—^-^Fi/\m'f,ml,mi,mi) 



+ 



>-4) 



4,0 ' 



m^,m^,mr^,m^) 



+ FiJ {ml, ml, ml, ml) + C(n - 4) 



j^(-2)/ 2 2 2 2n 

^4,0 {mi,m2,m^,m^) 



T7(-l)/ 2 2 2 2n 

^4,0 {^ii'n72,m^,m^) 



FlA{m^,m^,m^,m^ 



+ - 



V^^^ {ml, ml, ml) — V^^^ {ml, ml, ml 

2 2 

mf — ml 



The expansion in {n — 4) of the coefficient of p"^ in Eq. (jl6|) is 



(17) 



(1^ 



(19) 



F4 i{n, ml, ml, ml, ml) = C^{n)< 



T-i( — 2)/ 2 2 2 2\ 

Fli {mi,m2,m^,m^} 



{n - 4)^ 



1 



{n-A) 
with 

T7(-2)/ 2 2 2 2n 

Fl^ '{m^,m^,m,^,m^) 

T7(-l)/ 2 2 2 2\ 

^4,1 {mi,m2,m^,m^) 

17(0)/ 2 2 2 2\ 

Fl{{m^,m2,mr^,m^) 
ml{ml — ml — ml) 



^2^ ^3, ^l) +^4°i ("^1' "^25 ^1, , ^1) + 0{n - 4) 



Sf\ml,ml) 



ml 



{ml — m4)^i?^(mi, m^, m"^ 



7 + - log(mi) log 



ml 



+ ^og{ml) [log{m 



+ log(m^)^ (-3 + ^log(m^; 



^og^i^lml) - 31og(m2m3) 



(20) 



+ T log(mo 



6 



+ 



2 2 

mfm2 



/I 3 1 

7 + log(mi) f - - log(m2) + log(m3) "2 + 4 

+ log(m^) Q log(m2) - ^ + ^ log(m^; 



+ 



(m| — — mg) 
16 



7 + - log {m^m^) — 3 log(m]^m3) 



+ (rrig — + mj) V^*-^'' (m^, mg, mg) 



(m? - ml) [ 32 V 4 

,2 _2 _2 



- - + log(m4; 

(0) 



222 

, "^117(0)/ 2 2 2n I "^277(0)/ 2 2 2n I "^3 17(0)/ 2 2 



+ 



+ 



mi 



[ml — miy 4 1 8 



1(3 3 1 r 



mi 



ml)^ 



(m^, m^, ml) - V^(o) (m^, m^, m. 



(21) 



The above expression for F|°''(m^, m^, mg, m^), although not at a glance, is symmetric in 
the exchange of m| and mg. 

The coefficient of the double expansion, at n ~ 4 and ~ 0, Fq'i (m^, m2, mg) is explic- 
itly given in Eq.(51) of |3], while the other coefficients F^^^i (m^, m.|), F2°''(mf , m^, m|), 
F3°''(mf , ml, m|), can be easily found by using the relations Fj = —-^Fq, {i = 1, 2, 3). 



4 Numerics. 

As already illustrated for the case of the 2- loops sunrise graph in [TH], one can use the 
Runge-Kutta method fH] to advance the solution of the MDE, in this particular case 
Eg . (1121) . from the known initial conditions at ~ 0, following a complex p^-path in the 
lower half-plane. We add Eq. (jl2j) and Eq. (jl3|) to the system of the four MDE for the four 
sunrise Mi's jTH], as they are present also in Eq.lfT^ , and we solve the system at the 
same time for all the Mi's in the same numerical program. 

All the features discussed in fH] regarding precision and organization of the program 
remain unchanged. Of course the execution time increases as the system of equations is 
bigger. 

Again for convenience, we use reduced masses and a reduced external momentum 
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squared 

nii ^ 

mi^r = — , , Pr = ^, , /i = mi+m2 + m3, (22) 
/i /i"^ 

where the choice for /i is motivated by faster numerical calculations. 

As discussed in jTH], close to a special point, different fromp^ = 0, numerical problems 
arise. For F^'^\ml,ml,ml,ml,pl) the special points are: pf. = oo (analytically discussed 
in ini), Pr = 0, the 2-body threshold pf. = plbthr- — ~("^i,r +"^4,r)^, the 3-body threshold 
pI = pli,,th,r = — ("^i,r + rn2,r + '^s.r)^ and the 2-body and 3-body pseudo-thresholds 

Pr = Plb,ps,r = -("^l,r " "^4,r)^ pI = plb,psl,r = -("^l,r + "^2,r " ^«3,r)^ Pr = Plb,ps2,r = 
— {^l,r — ^2,r+^3,rY,Pl = P3b,ps3,r — ''{^l,r ~ TT^2,r — ^3,rY ■ To obtaiu ValueS of the MI's 

for p'^ close to thresholds and pseudo-thresholds we use the method discussed in detail in 
the next section. 

To illustrate some results of the program, and in the same time to show the behaviour 
of the function F4 over a wide range of the variable p^, we present here plots for two 
different sets of masses. In Fig. |2]we plot the values of Re F^^^ and Re S^^^ as functions 
of p,^, for the choice of the mass-values mi = 2, m2 = 1, 1713 = 4, = 20.1 (arbitrary 
units). The 2-body threshold, = —(mi + 7/14)^) , corresponds to pl^thr — —9.96755102. 
The presence of the 3-body threshold {plbthr = due to the very definition of the 
scale n) is hardly visible from the Re F^^ plot, but clearly from its imaginary part plotted 
in Fig|21 Its behaviour around the 2-body threshold is more 'fuzzy' due to the presence 
of root terms in the threshold expansion (see next section for an explicit form of the 
threshold expansion). The strikingly different behaviour of F^'^ and S^^^ near the 2-body 
threshold is due to additional logarithmic terms in the expansion of F^^ , which are absent 
in 5'*^°^ expansion. 

In Fig. m are plotted the values of Re F^^"^ and Re S^^^ as a function of pi, for the 
second choice of the mass-values, where one of the masses is considerably bigger than the 
others mi = 1, m2 = 9, = 200, = 20.1. This time the higher threshold is the 
2-body threshold at j-i-i^j, — —0.010095465, while the 3-body threshold is as always at 
Psbthr — Both functions show pronounced peaks around the 2-body threshold, and 
the 'zoom' of the peaks in Fig. shows that the qualitative behaviour of both functions 
nearby the 2-body threshold does not depend on masses (compare Fig. 

The imaginary part of Im F^^^ starts to be non vanishing at the 2-body threshold, as 
shown in Fig. El but this time no visible modification comes in at the 3-body threshold. 
The peaking behaviour of the imaginary parts is shown enlarged in Fig. [7| 

In Table ^ are reported the benchmark values of lie F^^\m1,m2, ml, ml, p^.) and 
Im Fl^\m1,m2,ml, ml, pI) for the masses mi = 1, m2 = 9, m^, = 200, m4 = 20.1, 
/X = mi + m2 + ma. 
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Pr 


Re fP 


Im 




-30. 


-0.44092793813(5) 


-0.38633149701(4) 




-15. 


-0.47753230147(5) 


-0.13848347339(5) 




-1.5 


-0.19215624308(2) 


0.5235669521(1) 


2 

P3b,th,r — 


-1.0 


-0.0990578391(1) 


0.60473101562(7) 




-0.99 


-0.09669632676(2) 


0.6066658942(1) 


P3b,ps3,r 


-0.981043 


-0.094558472357(8) 


0.60841511708(6) 




-0.9 


-0.07412586818(2) 


0.6249850792(1) 


2 

P3b,ps2,r~ 


-0.835918 


-0.056354400576(6) 


0.63914196276(5) 




-0.825 


-0.05316178728(1) 


0.6416578704(1) 


2 

P3b,psl,r~ 


-0.818594 


-0.051264529466(5) 


0.64314892655(5) 




-0.8 


-0.04565273027(1) 


0.6475413524(1) 




-0.1 


0.65325134562(1) 


0.99590695379(1) 


P2b,th,r = 


-0.010095 


2.101104208(4) 


0.0 


P2b,ps,r 


-0.008272 


1.697672304903(6) 


0.0 




0.0 


1.351553692518424 


0.0 




1.0 


0.21346230166(6) 


0.0 




30.0 


0.117099162908(6) 


0.0 



Table 1: The benchmark values of Re F4 (mf , m^, m|, m4,p^) and 

Im F^'^\ml,m2,'rn'l,m'l,p'^) for the masses mi = 1, m2 = 9, = 200, = 20.1, 

^ = mi + 1712 + ^3- 



The values at thresholds and pseudo-thresholds are obtained with the method dis- 
cussed in the next section, the value at = is from the analytic formula. 

5 Thresholds and pseudo-thresholds. 

The numerical calculation with the system of MDE does not allow for pi much closer than 
10""^ to the thresholds and pseudo-thresholds {special points)^ due to the numerical insta- 
bihty of the equations in these surrounding-domains. So a special treatment is required 
to get precise values of the Mi's there. 

For the sunrise Mi's the analytical expansions at pseudo-thresholds [22] and at thresh- 
old 121] were previously obtained. In [TH] these results were used as starting points for the 
system of MDE's for the sunrise Mi's to get reliable results in the surrounding- domain 
of these points. However this method is not universal, as it requires a separate analytic 
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calculation of the MI expansion at these points, which is difficult and probably not always 
possible. 

We discuss here in detail the features, advantages and limitations of an approximation 
method, which is rather precise, universal and easy to use. This method was recently 
proposed in 20 . However the 'universal' approximant suggested by the author applies 
only to some cases and fails at the 2-body threshold relevant for the present calculations. 

The method consists in the construction of a suitable approximant of the MI, which due 
to the smallness of the surrounding- domain is naturally the expansion of the MI around 
the considered special point. The proper form of the expansion (with undetermined coef- 
ficients) around each special point can easily be deduced from the MDE's themselves • 

The coefficients of the approximant are calculated using numerical values of the MI, 
obtained solving the system of MDE for some points nearby the special point, but outside 
its surrounding- domain to avoid numerical instability. The approximant is used to get 
the values of the MI at the special point and within its surrounding- domain. Note that 
the proper form of the expansion around a special point is crucial to produce right results 
around special points. With the wrong choice of the approximant one may still get the 
values at the special points right, when the points chosen for the approximant 'construc- 
tion' are symmetrically distributed around the special point. The values obtained for the 
surrounding- domains will be however wrong. 

To test the precision of this approximation method, we compare its result at the 
special point with that obtained from analytic result there, when available, and in the 
surrounding- domain with the result given by the advanced numerical solution of the 
MDE's started from the special point. Of course the approximation method can be used 
also in a region where a MI does not contain a special point to test the precision of the 
method. 

The 2 and 3-body pseudo-thresholds are the simplest points to be treated in this 
way: as they are regular points the approximant can simply be the power expansion. For 
each MI {fI^\ i = 0, 1,2,3,4 and S^^^) in each of its pseudo-thresholds can be used an 
approximant of the type 

F{x) = ao + aix + a2x'^ + a^x^ , (23) 

where x = p"^ — Pp^ ,^ and Pp^ ^. is the value of p"^ at the pseudo-threshold. As we want to 
use X well below lO"'^, in Eq.()23p it is enough to use four terms in the expansion, the 
truncation causing a relative error of the order of x^ = 10~^^. The actual implementation 
of F{x) uses the four points x = ±10~^, ±0.5 ■ 10~^ to calculate the four coefficients 
Qi, i = 0,1, 2, 3, the error associated to its value is dominated by the error coming from 
the numerical calculation of the MI at these points. It has to be noted that the coefficients 
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of the expansion Eq. can be complex if the expanded function is complex valued - as 
it is the case above threshold(s). 

The values of the sunrise Mi's Fj^^\ i = 0,1,2,3 at the 3-body pseudo-thresholds 
and of the 1-loop self-mass MI S^^^ at the 2-body pseudo-threshold, can be obtained 
from known analytical formulas [2S1 1^ • The approximation method can reach in the 
implemented program a relative error of the order of 10~^^ — 10~^^ and its results are 
always in agreement with the analytical ones within this error. 

In the surrounding- domain (where x < 0.5 • 10^^) of the 3-body pseudo-threshold 
approximation method results for the sunrise Mi's are in agreement with the results 
obtained from MDE advancing solution starting from the pseudo-threshold ^H] within a 
relative accuracy of the order of 10~^, which is the best accuracy obtainable by the MDE 
advancing solution program in this region. However, as the approximation method values 
for F^^\ i = 0,1, 2, 3, and 5'*^°'' reproduce well the values from analytical formulas at the 
pseudo-thresholds, and from MDE advancing solution program starting from = in 
all tested points outside the surrounding- domain (where x > 0.5 ■ 10"^), we are confident 
that the approximation method is better and that the stated relative accuracy of order 
10^^^ — 10~^^ is a safe estimate. 

At the 3-body threshold the proper approximants for the -F/"'', i = 0,1, 2, 3, 4 are 
Fj;°\x) = aifl + (ai,i + 6^,1 log(x)) x + (0^,2 + &i,2 log(x)) x^ + (0^,3 + 6^,3 log(x)) x^ , (24) 

where X = pi — plh^th,r pli,^th,r = ^("^l,r + "^2,r + "T-S.r)^ = " 1- Wc USC X = 

±0.5 ■ 10~^, ±10~^, ±1.5 • 10~^, ±2 ■ 10~^ to obtain numerically the coefficients ajj, i = 
0, 1, 2, 3, 4, j = 0, 1, 2, 3 and bi^k, k = 1,2, 3. Again all the coefficients in the expansion of 
^4°"* can be complex depending on the values of the 2-body and 3-body thresholds. For 
S^'^^ we use Eq.d^Hl as it is regular at the 3-body threshold. 

Comparisons with the analytical results for the sunrise MI at the 3-body threshold jTS] 
and other tests like those at pseudo-thresholds were performed with the same conclusions, 
so again a precision of the order of 10^^^ — 10^^^ can be assumed. 

At the 2-body threshold the proper approximants for S^^^ and F^'^^ are 
S^'^^ = ^5,0 ± o-saX + as,2x'^ + as,3X^ + \fx {bs,o + bs,ix + hs,2X^ ± hs^^x 



FX = 04,0 ± aA,ix + a4^2X + 04,3^ ± vx 
/x log(x 



C4,o ± C4,iX ± C4,2X^ ± 04,3^^ 



(25) 



where x = p^. — P2b.thr vibthr — ~{^i,r ± "^4,r-)^- The threshold expansion can be 
easily derived from Eq. ()12j) for fI^\ while for S^^^ it was given in |1^. The real constants 
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o-s^ji bsj and real or complex constants 04^, 64^,04^, j = 0, 1,2,3 (see discussion above) 
are obtained calculating numerically with the MDE advancing solution the values of S^^^ 
at X = ±^■10-^ X = ±i■0.5■10-^ x = ±i•0.25■10-^ x = ±z-0.125- 10"^, and of there 
and also at x = ±z- 0.0625- 10^^, x = ±z- 0.03125 ■10"'^. The choice of imaginary values for 
X allows to overcome the numerical problems arising in the present case for real values of 
x when |x| < 10~^ with the usual MDE advancing solution program starting fromp^ = 0, 
due to the steep behaviour of the functions S^^^ and F^^^ nearby the 2-body threshold, as 
can be seen in Fig. On the other hand the MI values at |x| < 10~^ are needed to get 
enough precise values for the coefficients in the approximants. For -P/^"*, [i = 0, 1, 2, 3) we 
use Eq.(I2ni) as approximants. 

For the S^^^ approximant values a relative precision of the order of 10^^^ — 10^^^ is 
obtained as before, confirmed by the comparisons with the analytic result. 

For the F^^^ approximant values the contrasting requirements does not allow a relative 
precision better than 10"^ — 10~^. 

6 Comparisons. 

Numerical values for the MI of Eq.(0) for non- vanishing values of the masses are reported 
in literature, only for small values of time-like p^, in [7j and more extensively in jS]. More 
precisely there was calculated the quantity 

Ti234Nip'^; rhj, rhj, m^, m^) = Ti234(p^ rhl, rhj, ml, rhl, f/) - Ti234(p^ rhl, rhl, 0, 0, /i^) 
= 16Fi°)K,m2,m2,m2,p^/i2) _T(°l4(-p2;^l^?,0,0,^2) + ^C(2) , (26) 

with T— functions defined in [HIH]. The function Ti234Ar(p^; m^, m|, m|, 7714) is finite in 
the n — >■ 4 limit and it is independent from /i^, the arbitrary mass scale introduced in 
Eq.Q. Here the /i^ dependence is explicitly shown, whenever present. The quantity 
Ti23a{p^', ml, ffi\, ffi\, fh\, /i^) is like the MI of Eq.(^, with a different normalization factor, 
which results in the factor 16 and the additive quantity C,{2)/2, the second line of Eq. pH|) . 
The analytical expression for Ti234(p^; fn{, fh^, 0, 0, /i^) is given in [2S]. 

Due to the different assignment of the mass indexes (m^) to the lines in our Fig^ 
with respect to the correspondent ones (m^) in [71 |H1 I^EI, we have the correspondence 
mi = Th2 = 3, m2 = m^ = 7, m^ = m^ = 5, m^ = fhi = 1, where the numerical values 
are for the comparison with Table 1 of jS]. One has also p"^ = —p^, due to Minkowski- 
Euclidean space transformation. In the second line of Eq. lf^ Ff\m\, m\, m\, m\,p'^, /x^) 
is calculated through the differential equation Eq. lfT^ and T^2^^{—p^;m\,m\,Q,Q,jjL^) is 
the finite part in the n — >• 4 limit of Ti234(— p^; m\, m\, 0, 0, /i^) from j^Hl- 
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^2 
P 


A 

A 


TD 

D 


-Z.bODD ( 


-8.45038098/50(01 j 


-8.45038 


-1. ( ( ( io 


-o.zo / 4 / ioiooz(^oz j 


-o.zo (4 / 


-1.18519 


-8.18480592994(33) 


-8.18481 


-0.790124 


-8.11877664655(31) 


-8.11878 


-0.526749 


-8.07577187774(28) 


-8.07577 


-0.351166 


-8.04753658726(33) 


-8.04754 


-0.234111 


-8.02890161217(30) 


-8.02890 


-0.156074 


-8.01656069505(26) 


-8.01656 


-0.104049 


-8.00836962582(21) 


-8.00837 


-0.069366 


-8.00292496915(41) 


-8.00292 


-0.046244 


-7.99930227781(36) 


-7.99930 


-0.030829 


-7.99689023488(21) 


-7.99689 


-0.020553 


-7.99528370180(10) 


-7.99528 


-0.013702 


-7.99421324511(10) 


-7.99421 


-0.009135 


-7.99349993358(10) 


-7.99350 


-0.006090 


-7.99302446227(10) 


-7.99302 



Table 2: Values of Ti234Ar(p^; m^? ""^s; '"^4) small values of = —p^. The masses 
are mi = ffi2 = 3, m2 = fhi = 7, 1713 = rh^ = 5, = rhi = 1. In column A are reported 
the results of the present method, with in brackets the error in the last digits, in column B 
are reported the results from column D of Table 1 of [Hj, to whom barred notation refers. 
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Our results are shown in column A of Table |21 and are in complete agreement with the 
results in column D of Table 1 of [Hj, obtained from a one-dimensional integral represen- 
tation and reported also here in column B for the reader's convenience. The results are 
purely real for small values of i.e. smaller than the smallest threshold value, which 
for the chosen values of the masses is the 3-body threshold of Ti2u{—p'^', ml, ml, 0, 0, /i^), 
corresponding to p"^ = —(mi + + 0)^ = —9. 

Recently in (201 ^ calculation with the same method has provided the MDE and the 
expansion at = for all the 2-loop self-mass MI. Some numerical results are presented 
there for the function U corresponding to fI^\ that we have verified up to the 9 digits 
reported in [201 ■ The comparison is by no means simple, due to the different definition of 
the pole-term in the 4-denominator function in |201- The finite terms of U and F4 differ 
and the term proportional to {n — 4) of the Laurent-expansion of the 1-loop self-mass MI, 
that we take from [2Z|, is necessary to account for that difference. 

7 Conclusions. 

We presented the results of the analytical expansion around = and of the numerical 
calculation of the MI related to the general massive 2-loop 4-denominator self-mass dia- 
gram, using the 4-th order Runge-Kutta advancing solution method applied to the MDE 
satisfied by the MI. 

An approximation method, suggested in [20|, was properly developed here, allowing 
for precise numerical evaluation of the Mi's even in the surrounding- domains of the special 
points. 

The relative precision reached by the developed numerical program is of the order of 
10^^^-10~^^ for all points, but the surrounding- domain of the 2-body threshold, where 
the relative precision is 'only' of the order of 10^^-10~^. 

It is to be emphasized that the method requires only the knowledge of Master Dif- 
ferential Equations (MDE's) satisfied by the Mi's. The initial conditions follow from 
requirement of the regularity of MI at p^ = and MDE, while the form of the approxi- 
mants used in the surrounding- domains of the special points follows from the MDE only. 

Acknowledgments. Henryk Czyz and Agnieszka Grzelihska are grateful to the 
Bologna Section of INFN and to the Department of Physics of the Bologna University for 
support and kind hospitality. 
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Figure 2: Plots of Re F^^^ {mf, m^, ml, ml, pi) and Re S^^^ {m\, m\, pi) as a function of pi 
for mi = 2, m2 = 1, m^ = 4, m4 = 20.1 and fi = mi + m2 + m^. 
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Figure 3: Plots of Im F^\m\, m\, m\, m\, p^) and Im S'^^\m\, m\, p^) as a function of 
pI for mi = 2, m2 = 1, ms = 4, 7714 = 20.1 and = mi + m2 + ms. 
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Figure 4: Plots of Re F4 (mf , m^, mg, 7714, p^) and Re S^^^ (ml, ml, pf.) as a function of p.^ 
for mi = 1, 1712 = 9, ms = 200, = 20.1 and fi = mi + m2 + m^. 
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Figure 5: Enlargements for the peaks of the plots of Re F^\m\,m\,m\,m\,p1) and 
Re S'^'^\m\,m\,p1) as a function of p^ for mi = 1, m2 = 9, m^ = 200, m^ = 20.1 and 
/i = mi + m2 + mg. 
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Figure 6: Plots of Im (m^, m^, ml, ml, p'^.) and Im S^'^^mf, ml, p^) as a function of 
pI for mi = 1, m2 = 9, ma = 200, = 20.1 and /i = mi + m2 + ms. 
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Figure 7: Enlargements for the peaks of the plots of Im F4 \m\,m'^,m\,ml,p1) and 
Im S'^^^ml, ml, p^) as a function of pi for mi = 1, m2 = 9, ms = 200, m^ = 20.1 and 
/i = mi + m2 + ma. 
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